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A many-flavor electron gas (MFEG) is analyzed, such as could be found in a multi-valley semi- 
conductor or semimetal. Using the re-derived polarizability for the MFEG an exact expression for 
the total energy of a uniform MFEG in the many-flavor approximation is found; the interacting 

energy per particle is shown to be — 0.574447(-Ehflo "i 4 )?i with E^ being the Hartree energy, 
ao Bohr radius, and m* particle effective mass. The short characteristic length-scale of the MFEG 
motivates a local density approximation, allowing a gradient expansion in the energy density, and 
the expansion scheme is applied to electron-hole drops, finding a new form for the density profile 
and its surface scaling properties. 
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I. INTRODUCTION 

For some semiconductors, at low temperatures and 
high density, electrons and holes condense into electron- 
hole drops, which provide a good testing ground for un- 
derstanding effects of electron-electron interactions [J. 
Some of the semiconductors (and also semimetals) that 
electron- hole drops form in [3, y], such as Si, Ge, and di- 
amond have conduction band minima near the Brillouin 
zone boundary, for example Si has six degenerate val- 
leys, see Fig. [TJ a Ge-Si alloy has ten degenerate valleys, 
and Pbi-^-ySn^MnyTe has twelve valleys in the S band 
[Jj. When the material is strained, valley degeneracy re- 
duces aja jj_8|jr 
iiill 
ley degeneracy could be regarded as a control parameter. 
Because of this, as well as degeneracy being large in some 
semiconductors, valley degeneracy might be a good pa- 
rameter with which to formulate a theory of electron-hole 
drops. 

Previous theoretical analyses of electron-hole drops 
B HI US HI [H, |H Hi used an expansion of the energy 
density with parameters found from separate energy cal- 
culations [9(. An alternative approach is to assume that 
each valley contains a different type of fermion, denoted 
by an additional quantum number, which we shall call the 
flavor, the total number of flavors (valleys) is v. Further 
motivation to study flavors stems from the fact that in 
some previous studies of multiply degenerate systems the 
number of flavors has not been well defined, for example 
heavy fermions [2|| [22], [28[ , charged domain walls [HI , 



u .which can be experimentally probed 
[ El El El El El , meaning that val- 



a super-strong magnetic field 0, and spin instabilities 
[30, [31| . Cold atom systems in optical lattices [32|, [33, [34[ 
have a well defined number of flavors but weak interac- 
tions between particles. In electron-hole drops however 
the number of flavors is well defined and interactions are 



strong. 

The ground state energy and pair correlation function 
of a free many-flavor electron gas (MFEG) were exam- 
ined using a numerical self-consistent approach for the 
local field correction by Gold [35|], and superconductiv- 



ity was studied by Cohen 
Keldysh and Onishchenko 



. Following the method of 
Andryushin et al. Q stud- 
ied the behavior of the free MFEG by summing over all 
orders of Green's function contributions, they found an 
exact expression for the correlation energy of a MFEG 
(which dominates the interacting energy in the extreme 
many-flavor limit). This paper describes the derivation 
of a more versatile formalism, based on a path integral, 
which gives an exact expression for the total energy of the 
MFEG; the theory could apply with as few as six flavors 
where the exchange energy assumed small by Andryushin 
et al. [3| would be significant. 

As well as studying the uniform case, the previously 
unstudied density response of a MFEG not constrained 
to be uniform is investigated. The screening length-scales 
of the MFEG are shown to be short relative to the in- 
verse Fermi momentum, suggesting that a local density 
approximation (LDA) might be a good approximation, 
motivating a gradient approximation. This gradient ex- 
pansion is then applied to analyze the electron-hole drop 
density profile, and to simulate effects of strain the scal- 
ing of drop surface thickness and tension with number of 
flavors is examined. 

In a MFEG with v flavors at low temperatures the 
relationship between the number density of electrons n 
and Fermi momentum pp is 



3tt 2 ' 



(1) 
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When the electrons have multiple flavors, each Fermi sur- 
face encloses fewer states so p? oc z; -1 ' 3 . The local band 
curvature governs the electron effective mass, the band 
structure is often such that the holes relax into a single 
valence band minimum at the T point (see Fig. [1]); here 
holes are assumed to be heavy and spread out uniformly 
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FIG. 1: The Si band structure in the [100] direction generated 
with a LDA-DFT approximation by a plane-wave pseudopo- 
tential method [33] . The Fermi energy is as E = eV; below 
are valence bands with holes at H; above conduction bands, 
the bold parabolic curve signifies the first conduction band 
valley with electrons at E. 



providing a jellium background. 

The Thomas- Fermi approximation predicts a screening 
length k^ 1 = (47re 2 g)~ 1 / 2 , where g is the density of states 
(DOS) at the Fermi surface. The DOS is dependent on 
the number of flavors as g oc v\JEp oc is 2 ' 3 and so k _1 oc 
j / -i/3_ rpjjg ra ti of the inverse Fermi momentum length- 
scale to the screening length varies with number of flavors 
as Pf/k oc v~ 2 ' 3 . This paper takes the many flavor limit 
v ~^> 1, in which the screening length is smaller than 
the inverse Fermi momentum, kT 1 < p F , the many- 
flavor limit therefore means that the wave vectors of the 
strongest electron-electron interactions obey q 3> pp, this 
is the opposite limit to the RPA which assumes that pp 3> 
k. Physically this means the characteristic length-scales 
of the MFEG are short so a LDA can be used in Sec. IIIII 
to develop a gradient expansion. 

The conduction band energy spectrum is characterized 
by two spectra E(q) and €j(p) as shown in Fig. [2] There 
are two energy functions: E(q L ) gives the energy in the 
band structure at momentum q; Cj(p) = p 2 /2m denotes 
the kinetic energy at momentum p with respect to the 
center of the ith valley (the dispersions of all valleys are 
assumed to be the same and isotropic so that €j(p) = 
e(p). Andryushin et al. [3| have outlined a method of 
calculating a scalar effective mass for anisotropic valleys) . 

Physical manifestations of the many-flavor limit in- 
clude effects where the DOS at the Fermi surface (energy 
Ep) is important; from Eqn. (fT]) the DOS of a particu- 
lar flavor i shrinks as gi(Ep) ex pp ex i/ -1 ' 3 whereas the 
DOS of all flavors grows since g(Ep) = Y^i=i9i(^F) = 
vgi(Ep) oc v 2 ! 3 . With increasing flavors, more electrons 




FIG. 2: The dark grey ellipsoids show Fermi surfaces of elec- 
trons in the six degenerate conduction band valleys in Si. 
■E(q) is the energy with momentum q measured with respect 
to the r point, e;(p) is energy with momentum p measured 
with respect to the center of the ith valley. 



are within ~ kp,T of the Fermi surface hence are able to 
be thermally excited, therefore the heat capacity of the 
MFEG, C = 12klT(v/3TT 2 n) 2 / 3 /5, increases with num- 
ber of flavors. The Stoner criterion [38|, [39[ for band fer- 
romagnctism states that for opposite spin electrons inter- 
acting with positive exchange energy U, ferromagnetism 
occurs when g{Ep)U > 1. With increasing number of 
flavors the total DOS g(Ep) oc v 2 ! 3 increases so that 
the Stoner criterion becomes more favorable. However, 
this analysis does not take into account the curvature of 
the DOS at the Fermi surface which can be an impor- 
tant factor in determining whether ferromagnetism oc- 
curs j40l . |41| . The effect of the total DOS is also seen 
in the paramagnetic susceptibility, this is proportional to 
the total DOS at the Fermi surface so is expected to in- 
crease with number of flavors. Analogously one can com- 
pare a transition metal with narrow d-bands that leads 
to a large DOS at the Fermi surface with a simple metal 
that has broader free electron conduction bands and so a 
lower DOS at the Fermi surface. Similar to many flavor 
systems, transition metals are experimentally observed 
[42 ] to have a significantly higher specific heat capacity 
and greater magnetic susceptibility than typical simple 
metals. The simple scaling relationships with number of 
flavors for heat capacity and magnetization provide ad- 
ditional motivation to analyze a MFEG in more detail. 

This paper uses the atomic system of units, that is 
e 2 = h = m = l/(47reo) = 1, but is modified so that 
m denotes an appropriate effective mass for the electron- 
hole bands, which is the same for all valleys. This mass 



m = m c m* can be expressed as a multiple of the elec- 
tron mass m c and the dimensionless effective mass m*. 
The units of length are then ag = ao/m* where ao is 
the Bohr radius, units of energy are those of an exciton, 
_Ej* = E^m* , where E^ is the Hartree energy. These six 
quantities, defined to be unity, give the standard atomic 
units when m* = 1 . For the important relationships that 
are derived in this paper particular to the MFEG, the 
full units are shown explicitly for clarity. Throughout 
this paper density is denoted by both n (number density 
of particles) and r s (Wigner-Seitz radius). 

In this paper, firstly a new formalism for the uniform 
system is derived. In Sec. II Al the system polarizability is 
found, in Sec. Ill Al the general quantum partition func- 
tion is derived, and in Scc. lIIBl thc uniform MFEG total 
energy is calculated. Secondly, we examine the system 
with non-uniform density: in Sec. IIIII a gradient expan- 
sion in the density for the total energy is found and is 
applied to electron-hole drops in Sec. HVi whose density 
profile and surface properties are calculated. 



A. Polarizability 

In this section the MFEG polarizability is derived. 
Though the result for the polarizability is the same as 
previous work [3, |43[, the derivation is presented here 
since an assumption made leads to an applicability con- 
straint on the many-flavor theory (in Sec. IIICj) . and the 
MFEG polarizability is an important quantity that will 
feature prominently in two main results of this paper: the 
Sec. |TT] derivation of the MFEG interacting energy and 
the derivation of the gradient expansion (see Sec. IIII|) . 

The polarizability IIq IF (q, u>) , where the superscript 
"MF" (many-flavor) denotes this is only for a MFEG, 
at wave vector q and Matsubara frequency to is given by 
the standard Lindhard form 



results in 



nJT 



■^— ' 47T 3 



St j f n F (ei (p)) - n F (ej (p + q )) 
e, (p) - €j (p + q) 



(qjW ) = y ^l / »*;^^-"^^T^ dP) 



(2) 

where n-p(ti) = l/(e^( e * ^ + 1) is the Fermi-Dirac distri- 
bution, j3 = 1/kftT, and fi is the chemical potential. In 
the standard expression for the polarizability the Fermi- 
Dirac distribution would contain the energy spectrum 
E(p) but in the MFEG each electron is in a particu- 
lar valley so the polarizability should be re-expressed in 
terms of the energy dispersion of each valley Cj(p) (see 
Fig. ^ and the contributions must be summed over the 
valleys i and j. Eqn. ([13]) shows that a large Coulomb 
potential energy penalty V(q) ex 1/q 2 inhibits exchange 
between different valleys so that the Kronecker delta Sij 
removes cross-flavor terms, and as all of the conduction 
valleys have the same dispersion a factor of v will replace 
the remaining summation over valleys. Supposing each 
conduction valley has a locally quadratic isotropic disper- 
sion relationship (with effective mass m), symmetrizing 
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dp. 

(3) 

In Sec. Ill Bl it is shown that the typical momentum ex- 
change q ~ (ha Q )n 1 / 4 is large relative to the Fermi 
momentum, therefore the two volumes in momentum 
space of integration variable p, defined by the two Fermi- 
Dirac distributions are far apart relative to their radii 
q/2 ^> p F , and the temperature is sufficiently low so that 
the high energy tails of the two distributions have neg- 
ligible overlap. Within these approximations the simple 
form for the polarizability is 



nr(q,«) = - 



(w/g) 2 + g 2 /4' 



(4) 



this expression agrees with the many-flavor polarizability 
found by Andryushin et al. |3| and Beni and Rice [431 ]. 

The standard Lindhard form for the polarizability, 
when taken in the same q 3> pf limit as imposed by the 
many-flavor system agrees with Eqn. (H|). In the static 
limit where frequencies are small compared to the mo- 
mentum transfer q 2 3> (H 3 / a^ 2 E£)ui , the polarizability 
varies as q~ 2 . In this limit, one Green's function is re- 
stricted by the sum over Matsubara frequencies to lie 
inside the Fermi surface whilst the other gives the polar- 
izability dependence of l/e(q) ~ q~ 2 due to the excited 
electron's kinetic energy. 

The derivation of the polarizability accounted only for 
intra-valley scattering, which means that in the MFEG 
the same terms contribute Q as in the RPA for the stan- 
dard electron gas. Therefore, diagrammatically, in the 
polarizability all electron loops are empty, the polariz- 
ability contains only reducible diagrams, which is de- 
noted by the polarizability subscript "0" . 



II. ANALYTIC FORMULATION 

Having reviewed the derivation of the polarizability of 
the system it is now used to formulate two complemen- 
tary components of the many-flavor theory. The first 
is the derivation of the energy of a uniform system; we 
begin by calculating the general quantum partition func- 
tion in Sec. Ill Al and continue for the homogeneous case 
in Sec. Ill Bl The validity of the many-flavor approach 
for a uniform MFEG is investigated in Sec. Ill CI The 
second part of the formalism is a gradient expansion of 
the energy density, looked at in Sec. IIIII Finally, the 
uniform and gradient expansion parts of the formalism 
will be brought together to study the model system of 
electron- hole drops in Scc. lIVI 



A. Partition function 

To derive an expression for the total energy of the sys- 
tem a functional path integral method is followed, which 
is a flexible approach that should be extendable to inves- 
tigate further possibilities such as modulated states and 
inter-valley scattering. Fcrmion field variables ip are used 
to describe the electrons irrespective of flavor in the dis- 
persion E(p). Overall the system is electrically neutral, 
so in momentum representation the q = element is 



ignored. The repulsive charge-charge interaction acting 
between electrons is V(r) = e 2 /r, we explicitly include 
the dependence on electron charge e (even though it is 
defined to be unity) so that the charge can be set equal to 
zero to recover the non-interacting theory. For general- 
ity we consider stationary charges Q(r) embedded in the 
MFEG which have a corresponding static potential f (r). 
The quantum partition function for the MFEG written 
as a Feynman path integral is then 



J 



Z = 



cxp ( / / V(r, t) (-iu + E (p) - n) i'{r, r)drdr J 



[ 2 



(ip(r', r)V>(r', r) - Q(r')) V(r - r') (ip(r, T )ip(r, t) - Q(r)) dr'drdr ) T>ipVtp. 



(5) 



This expression for the quantum partition function differs from that used for an electron gas (which has just a single 
flavor) only by the operator E(p) which gives the appropriate energy dispersion. To recover the standard electron 
gas result, which has a free particle dispersion relationship centered at the T point, one should set E(p) = p 2 /2m c . 
For the MFEG, as outlined in Fig. [21 E(p) represents the dispersion relationship of the whole conduction band, but 
no approximation concerning the flavors has yet been made, so the formalism applies for any number of flavors with 
a suitable energy dispersion relationship. 

To make the action quadratic in the fcrmion variable ip, the Hubbard-Stratonovich transformation [44| introduces 
an auxiliary boson field <fi(r, t) 



t 



/cxp I — Y^ <t>{l,u)(q 2 /4:ir)<t>(-(i, 



q#0,w 



, ire 2 8n v-^ 4ir 



cxp 



ip(r, r) f-iu) + E (p) + U - fi + ie$\ tp(r, r)drdr J VipVipVcp. 



(6) 



The direct decoupling channel [44| was chosen as the rel- 
evant contributions come from a RPA-typc contraction 
of operators. 

The term labeled with a (f ) exclusive of both fermion 
variables tp and the auxiliary field 4>, physically removes 
the electron self-interaction included when expressing 
the auxiliary field in a Fourier representation. Integrat- 
ing over the fermion variables tp and using ln(detA) = 
tr(lni) gives Z = J e^^^Vcp, where the action S[(p] is 



s\. 



pm 

~2^ 



J2 0(q,w)(q 2 /47r)^(-q,-w) 



q#o,w 



2ire 2 /3n ^ — 



q#0 



tx(lji(-iw + E(f>)+ie$ + U-fj,)). (7) 



g: 



Due to its similarity to an inverse Green function, GT is 
used to denote the argument of the logarithm, the sub- 
scripts "0" or "0" denote whether the inverse Green's 
function includes the auxiliary field or is free. 

Finally, we note for use later that the ground state 
total energy per particle -Eg = ■E'int + E$ can be split into 
two components. The interacting energy is -E^t (found 
in Sec. Ill B|) and the non-interacting energy is 



3 

To 



2/3 



(8) 



which is the energy with interaction between charges 
switched off (e = 0). It falls with increasing number 
of electron flavors due to the shrinking Fermi surface. 



B. Homogeneous Coulomb gas 

So far, up to Eqn. (J7J), the formalism is exact, however 
to perform the functional integral over bosonic variable <p 



an approximation must be made. To proceed one notes 
that with no external potential U(v) = the saddle- 
point auxiliary field of the action (Eqn. 0) is <f> = 0, 
fluctuations in the action are expanded about the saddle 
point solution in <f> giving the expression 



X 



SW\ 



ln(G^) + tr^F- 1 ^ - \ $G Q $Go 
i tr (^G^Go^Go^Go) + 0(0 6 ) 



e 2 4,,- 

q#0 H 



(9) 



Terms are now kept to quadratic order in <f>, this is equiv- 
alent to the RPA, analogous to the terms kept in the 
derivation of the many-flavor polarizability, see Sec. II Al 
This approximation will place a constraint on the validity 
of the formalism that is further examined in Sec. Ill CI 

The product of two Green's functions in the quadratic 
term in cf>, labeled (J), is identified with the polarizability 
Ilo , this is still expressed in terms of the general energy 
spectrum E(p) so is not yet necessarily many- flavor and 
does not carry the superscript "MF" used in Sec. II Al 
Following a multi-dimensional Gaussian integral over the 
fluctuating field </> (to quadratic order) the quantum par- 
tition function is 



Y[(q 2 /4n-e 2 Ilo(q,u;)y 



v^ 47re 2 

1 ^ £— " 



(10) 

In the low temperature limit we consider the free energy 
to get E G = -lim /3 ^ 00 (ln(.2)//?) = E + E int to get 
the interacting energy per particle, normalized so that 
Ei n t — with no interactions (e = 0), 



E„ 



1 



47re 2 \ dwdq 



e n 




(11) 



This equation remains general and is not necessarily in 
the many-flavor limit, it is in agreement with previous 
expressions for the interacting energy [5| studied not in 
the many-flavor limit, but which use alternative forms for 
the polarizability. If the standard (single flavor) electron 
gas form for the polarizability, the Lindhard function, is 
used then it is possible to recover, in the high density 
limit, the Gell-Mann Bruckner [45| expression for the to- 
tal energy 

However, to proceed, one should now assume many- 
flavors and use the appropriate polarizability, Eqn. ^. 
The many-flavor polarizability summed over all Matsub- 
ara frequencies in the zero temperature limit (3 — > oo sat- 
isfies i ~^2 U IIq (q, a;) = — n. This is used to substitute 



for the electron density n in the final term in Eqn. (fTTj) 
to yield the many-flavor result 
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(12) 



To evaluate the interacting energy one first substitutes 
for the many-flavor polarization using Eqn. ((][]), makes 
the change of variables fl = uj/q 2 and Q = q/n 1 / 4 , and 
re-arranges to get 



£int = -™ V4 
1 



167re 2 ,, / 16^e 2 /Q 4 \ 
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(13) 



The integral is independent of density and num- 
ber of flavors, so is the numerical factor ^4 3 rj = 
(^a; 3/4 )r(-5/4)r(3/4)/(2^ 5 / 4 ) « 0.574447(^a; 3/4 ) 
that was evaluated analytically [551 ]. The interacting en- 
ergy is therefore 



E- mt = -Awn 114 , 



(14) 



which is independent of the number of flavors. In eval- 
uating Eqn. (|13[) the main contribution to the integral 

over Q = q/n 1 ' 4 is at a momentum q oc {ha* Q )n 1 ' , 
so the interaction and screening length-scale in a MFEG 
is A ~ H/q ex a* n~ 1 ' 4 <C fr/pF, which is shorter than 
the Fermi momentum length-scale. 

The interacting energy E- m t = E cx + E COTI can be split 
into exchange energy E cx and correlation energy E coy:l . 
The interacting energy is independent of number of fla- 
vors, the exchange energy, E cyi = — (3/2)(3n/7ri/) 1/ ' 3 [3a |. 
falls with number of flavors, therefore the correlation en- 
ergy dominates over the exchange energy in the inter- 
acting energy in the many-flavor limit. In terms of the 
total energy, the correlation energy also dominates over 
the non-interacting energy, the kinetic energy that falls 
with number of flavors as Eq oc v~ 2 t z . The increas- 
ing importance of the correlation energy can be under- 
stood further by considering the electron pair correlation 
function. With increasing number of flavors the length- 
scales between electrons of the same flavor increase as 
ex iA> 3 r s and thus exchange energy and kinetic energy re- 
duce whereas the correlation energy depends only on the 
distance r s between electrons so is unaffected by the num- 
ber of flavors present. Andryushin et al. [3j and Keldysh 
and Onishchenko [2J found the Eqn. (fT4|) expression to be 
the correlation rather than interacting energy, neglecting 
the exchange energy which is small in the extreme many 
flavor limit. In Sec. Ill CI the expression for the inter- 
acting energ y is compared with self-consistent numerical 
calculations [35| on a MFEG with up to six flavors. 



The interac ting energy of a standard electron gas with 
a single flavor [35[ is more negative than that of a MFEG, 
which in turn is more negative than that of a Bose con- 
densate J46| , this could be due to the reducing negativity 
of the exchange energy, important in the single flavor 
system, but zero in the Bose condensate. In these two 
extreme systems, the single flavor electron gas and the 
Bose condensate, there is no notion of valley degeneracy, 
and therefore the intermediate system, the MFEG, might 
be expected at most to have only a weak dependence on 
number of valleys. In fact, the interacting energy of the 
MFEG, over the range of density found in Sec. Ill CI con- 
tains no dependence on the number of valleys. The ab- 
sence of flavor dependence is also present in the universal 
behavior for the exchange-correlation energy in electron- 
hole liquids proposed by Vashishta and Kalia J47I ]. 

The non- interacting energy term Eq oc (n/v) 2 ' 3 fa- 
vors low electron density, the interacting term E m t = 
— A^D-n 1 ^ favors high electron density, therefore the to- 
tal energy per particle has a minimum as a function of 
density of i?Gmin oc —v 2 ^ at n m ; n oc v 6 / b . The presence 
of a minimum in energy with density of the MFEG is 
consistent with the results of Andryushin et al. [3| and 
Brinkman and Rice [5[ who analyzed conduction elec- 
trons in a semiconductor. One consequence of this mini- 
mum is the possibility of a low density phase coexisting 
with excitons. 

Before analyzing the non-uniform system in detail in 
Sec. IIIII we can make qualitative arguments about its ex- 
pected behavior within a potential well. According to 
Thomas-Fermi theory, an electron gas in a slowly vary- 
ing attractive potential has a constant chemical potential. 
The electron gas is least dense at the edges of the poten- 
tial and is densest at the center of the well. In a MFEG, 
due to the negative interacting energy E- lnt = — A^n 1 ^ 
favoring high electron density, the density is expected to 
further reduce at the edges of the attractive potential and 
increase at the center of the well. In a repulsive potential 
the opposite should occur. 



C. Density limits 

In this section we will derive approximate expressions 
for the upper and lower density limits over which the 
many-flavor limit applies, these will be used to check the 
theory against numerical results [3a ] and to predict a 
lower bound on the number of flavors required for the 
theory to apply. 

To find the upper density limit one notes that Eqn. QI3p 
implies that an acceptable upper limit to the momen- 
tum integral would scale as q — a(fia )n l ' i , the 
constant a ss 4 was determined numerically and was 
chosen to give the q upper limit on the integral that 
recovered 95 % of the interacting energy. Addition- 
ally, the two regions of integration defined by the Fermi 
Dirac distributions in Eqn. ^ must not overlap, requir- 
ing that q/2 > pp. Combining these requires that for 



the many-flavor limit to apply the density must satisfy 
na* Q 3 < (a 12 i/ 4 )/(2 12 3 4 7r 8 ). Physically the breakdown at 
high density is due to the strongest interactions taking 
place on length-scales longer than the inverse length pp 1 . 

The low density limit is derived by considering the 
expansion of the action in the auxiliary boson field </>, 
Eqn. ©• In order to evaluate the Gaussian functional 
integral over the bosonic variable 4> it is necessary to ne- 
glect the quartic term in cf>, valid only when investigating 
the system with respect to its long-range behavior, that 
is pF^aS ^> 1 [3i an d therefore na^ 3 ^$> v/Sir 2 . The 
breakdown at low density can be understood because the 
MFEG is effectively a boson gas, all electrons will be in 
the T state (k = 0) and there is no exchange energy. 

The upper and lower critical density limits can be 
combined to conclude that the many-flavor limit result 
for interacting energy applies for densities that obey 
0.03^ <§C na,Q 3 <C 0.005^ 4 ; this density range increases 
as v 4 with number of flavors, the scaling relationship is 
the same as the applicable density range of the correla- 
tion energy found by Andryushin et al. [3j, though they 
did not provide estimates of numerical factors. 

Using the above high and low density limits it is possi- 
ble to estimate the minimum number of flavors required 
for the theory to apply. This is done by setting the lower 
and upper estimates for the allowable density to be equal, 
which gives v jj 2, this estimate is approximate due to the 
possible inaccuracies in the upper and lower critical den- 
sities used in its derivation. As the upper and lower crit- 
ical densities have been set equal, the many-flavor theory 
will apply here only over a very narrow range of densities, 
but this range widens with increasing number of flavors 
as i/ 4 . An alternative limit can be found by comparing 
the interaction energy predicted by the theory over the 
expected density range 1 of applicability with the results 
of Gold [3a]. Their numerical self-consistent approach 
gives interaction energies accurate to approximately 3 % 
when compared with single-flavor electron gas quantum 
Monte Carlo (QMC) calculation s I48l . [49J ] and some initial 
many-flavor QMC calculations (5CJ. At two flavors the 
interacting energy predicted by the many-flavor theory 
is ~ 10 % more positive than the self-consistent numer- 
ical results [351 ] indicating the many-flavor theory does 
not apply at two flavors. For six flavors over the pre- 
dicted allowed density range the many-flavor theory is 
between ~ % and ~ 4 % more positive than the numer- 
ical results, which indicates that the many-flavor theory 
can be applied within the predicted range of applicability 
(0.5 < r s /a,Q < I). The theory should be applicable in 
common multi-valley compounds, such as silicon which 
has six conduction band valleys, and to those with more 
valleys [J]. This result is corroborated by the results of 
initial QMC calculations [5(| on systems with between 6 
and 24 flavors. 

In the first half of this paper a new versatile formal- 
ism to describe a MFEG has been developed that could 
apply in systems containing approximately six or more 
degenerate conduction valleys. An exact expression for 



the total energy of the uniform MFEG was found and 
the applicable density range derived. The next step is to 
investigate the response of the MFEG to an external po- 
tential. A gradient approximation is developed in Sec. lIIII 
and this is applied to electron-hole drops in Sec. |lVl 



III. GRADIENT CORRECTION 



q/2 <C pp, reducing the volume available for integration 
so reducing the coefficient to 1/72. 

The gradient correction for the energy could be used in 
analytical approximations or as a DFT functional. This 
energy density expansion allows the MFEG to be applied 
to a variety of systems, its use for studying electron-hole 
drops is demonstrated in Sec. HVl 



In Sec. U it was shown that the typical length-scales of 
the MFEG are short q ^> p-p, motivating a local density 
approximation (LDA). This motivation is in addition to 
the usual reasons for the success of the LDA in density- 
functional theory (DFT) [5l[ - that the LDA exchange- 
correlation hole need only provide a good approximation 
for the spherical average of the exchange-correlation hole 
and obey the sum rule |52||. In this section the LDA is 
used with the polarizability derived in Sec. II Al to develop 
a gradient correction to the energy density that allows the 
theory to be applied to a non- uniform MFEG. 

The typical momentum transfer in the MFEG is q ~ 
(haQ~ )n 1 / 4 hence the shortest length-scale over which 

the LDA may be made is approximately (oq /K)n~ l ' A 
and the maximum permissible gradients in electron den- 
sity are |Vn| ma x ~ qn ~ (^ a o~ ))n 5 ^ 4 . The gradient 
expansion will break down for short scale phenomena, 
for example a Mott insulator transition. To derive an 
energy density gradient expansion we follow Hohcnbcrg 
and Kohn [53| and Rice [2l| and consider an external 
charge distribution n ex t(q) that couples to the induced 
charge distribution nj n d(q) with Coulomb energy density 



1 v-^ 47re 2 , . 
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(q)- 



(15) 



One now substitutes for n ext (q) using the relative permit- 
tivity l/e(q) = 1 + 7i ind (q)/n ext (q) = 1/(1 - 47rn^ F /g 2 ) 
and the many- flavor polarizability Eqn. (U|). The highest 
order term in l/q 2 gives the induced charge Coulomb en- 
ergy, the term of order q 2 is associated with a gradient 
expansion, in real space this gives the expansion for the 
total energy per particle 
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here Eq is ground state energy of the uniform system 
found in Sec. Ill Bl The form of the energy correction is 
similar to the von Weizsacker term [5J] , although here it 
is larger, having a coefficient of 1/8 rather than 1/72 as 
in the von Weizsacker case. The difference can be qual- 
itatively understood by considering the Fermi surfaces 
involved in the two cases for a given wave vector q, in 
the many-flavor case the Fermi surfaces involved in the 
integral of Eqn. $3§ do not overlap as q/2 ^> pf so there 
is a large volume in Fermi space available hence a large 
coefficient of 1/8, whereas in the ordinary electron gas 
(single flavor) the same Fermi surfaces do overlap, as now 



IV. ELECTRON-HOLE DROPS 

In this section the MFEG is applied to a simple sys- 
tem to investigate the properties of, electron-hole drops. 
An electron-hole drop is a two-phase system: a spheri- 
cal region of a MFEG surrounded by an exciton gas [43| . 
The density profile, surface thickness and surface tension 
of drops are investigated; the scaling of surface thickness 
and tension with number of flavors is also found since 
this is can be experimentally probed through externally 
imposed strain reducing the valley degeneracy [7J]. 

There have been four main theoretical methods used 
to analyze an electron-hole drop in silicon and germa- 
nium, semiconductors which have six and four flavors re- 
spectively. Rice [22| fitted an analytic form to the energy 
density minimum and included the lowest order of a local 
gradient correction, from the equation for energy density 
an analytic form for the density profile was derived. A 
similar approach was used by Sander et al. [201 ] and Rice 
j2l| to study the surface structure in more detail. The 
second approach [23|, [24[, which was also applicable to 
situations with an external magnetic field and uniaxial 
strain, conserved momentum, particle number and pres- 
sure balance at the drop surface, the resulting equations 
were then solved numerically. A third approach followed 
by Kalia and Vashishta [g| used a Pade approximant for 
the energy density [9|] derived specifically for silicon and 
germanium but did not include a gradient correction fac- 
tor. The fourth approach of Reinecke et al. [25[ again 
used a Pade approximant for the energy density and also 
included a gradient correction factor. The latter two ap- 
proaches assumed an exponential density profile for the 
drop. These four methods all use approximate forms for 
the energy density, an advantage of the many-flavor ap- 
proach is that the exact form for the analytic energy den- 
sity (within the many flavor assumption) can be used to 
solve for the drop density profile. Whilst analytic forms 
for the inner and outer density profile as well as a model 
for the entire profile can be derived, the general problem 
must be solved numerically. With an exact form for the 
density profile, electron-hole drop surface effects can be 
studied. 

Local charge neutrality is assumed so that the density 
of electrons and holes are everywhere identically equal. A 
LDA with gradient correction is used so the drop energy 
density is written as the sum of the local non-interacting, 
local interacting and the lowest order term in a gradient 



expansion, 



A 3D n(r)^ + 



5/4 | (V"(r)) 2 



8n(r) 



(17) 

The total energy of a drop is / e(r)dr and the total num- 
ber of electrons in the drop is J n(r)dr. The total energy 
is minimized with respect to electron density n(r) whilst 
keeping a constant number of electrons in the drop by 
applying the Euler-Lagrange equation with a Lagrange 
multiplier /i, which represents the chemical potential. If 
the drop has spherical symmetry the density must satisfy 
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rn 8/3 - 4(L4 3D rn 9/4 - 32/irn 2 {18) 



The boundary conditions are specified at the center of 
the drop, where the density takes the equilibrium homo- 
geneous MFEG value and the density is smooth, namely 
n(0) = n and n'(0) = 0. The differential equation 
Eqn. (fT8|) cannot be solved analytically for n(r), but a 
solution, n(r) = n, exists for fi = which corresponds 
to the homogeneous MFEG, that is a drop containing 
an infinite number of electrons. Before solving the dif- 
ferential equation numerically, two approximate schemes 
are developed, one that applies near the drop center and 
the other near the drop edge, and their predictions are 
compared with existing density profile forms. 

Near the center of the drop a perturbation solution 
about the equilibrium density, n(r) — h + An(r) where 
An(r) <C n, is considered. The solution to Eqn. (fT8"|) for 
the density is then 



n(r) = n+ — r- 1 
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sinh(Qr) 
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This density profile is characterized by an exponential 
reduction of the density away from n at the center. The 
energy Q 2 is physically the rate of change of energy per 
unit volume with respect to changing particle density, 
with 
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The second approximation scheme applies in the drop 
tail where electron density is low, n(r) <C n. The term 
containing the chemical potential is disregarded as it is 
arbitrarily small for the large drops under investigation, 
the non-interacting and interacting energy terms contain 
higher powers of density so are negligibly small. In this 
regime the solution to Eqn. (I18[) is 




r < r , 
r > r . 



(21) 



Here no and rg are variational parameters which must be 
fitted to a numerical solution. This analytic form shows 
that the electron- hole drop has a definite outer radius r , 
which is approached parabolically, it is also noted that in 
the drop tail the solution obeys the differential equation 
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If electron density is mapped onto a wave function tp 
through n = \ip\ 2 then the solution Eqn. (f2"Tj) obeys 
Schrodinger's equation at low energy, that is V 2 ?/> = 0. 
The implied Schrodingcr equation is for a low density 
MFEG with negligible interaction between electrons due 
to their large separation, consistent with the original as- 
sumption of low density in the drop tail. 

Previous studies of electron-hole drops [3, [22|, [25( had a 
solution with the same exponential form both inside and 
outside of the drop, our inner functional form, an expo- 
nential, agrees with previous work [8|, [22|, [25[, but our 
outer functional form, a quadratic-like polynomial, does 
not agree with the exponential decay seen in previous 
work. However, at the outside of the drop density is low 
and the arguments of Sec. Ill Cl show the many flavor the- 
ory, which requires that the density satisfies n 3> 0.03i/, 
does not apply here. The other theories [8|, [23, [25[ also 
do not apply in the low density region so both the many- 
flavor and previous theories fail to agree only where they 
are not applicable. 

Using just the solution for the density in the drop tail 
Eqn. (|21[) , a reasonable analytical approximation for the 
density form of the whole drop is 



n(r) = 



1 



1 



n (l/r - l/r ) 



(23) 



This solution has the correct functional form at both the 
inside (n(r) — > n) and outside of the drop and extrap- 
olates smoothly in between. It can be fitted to the ac- 
tual solution using parameters no and rg. However, the 
general differential equation is solved numerically giving 
the density profile shown in Fig. [3l The numerical solu- 
tion is well approximated in the inner and outer regions 
by Eqn. (TT9"|) and Eqn. (|21|) respectively, and the model 
Eqn. (|23|) provides a good fit to the numerical solution, 
having just a slightly too shallow gradient around the me- 
dian density but it agrees at both the center and outside 
of the well. 

To allow us to compare properties of electron-hole 
drops predicted using many-flavor theory with other work 
[a [25( , one can characterize the electron-hole drop prop- 
erties through its surface thickness D and tension 7. The 
surface thickness D is the width over which the density 
falls from 90 % to 10 % of its homogeneous equilibrium 
value n. The total surface energy is the difference be- 
tween the energy per unit area of the MFEG in the drop 
and the energy of the same number of particles at equi- 
librium density in a homogeneous system. The surface 
tension 7 is the total surface energy divided by the char- 
acteristic drop surface area, here taken to be the area of 
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FIG. 5: The variation of surface tension (dashed line and 
crosses) and surface thickness (solid line and pluses) with 
number of flavors present, the straight line fits are used to 
give the exponents of the scaling parameters. 



FIG. 3: The density profile of a 12 flavor electron-hole drop 
with density parameter r s — 1. The numerical solution is 
shown by the dotted line, analytical approximations to the 
inside (outside) of the drop by the dashed (dot-dashed) lines, 
and a best fit model fitted to the numerical solution by the 
solid line. 
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FIG. 4: The surface tension 7 of the 12 flavor drop of radius 
r m is shown using the solid line and pluses based on the left- 
hand y-axis. The variation of the surface thickness D is shown 
using the dashed line and crosses based on the right-hand axis, 
each point represents a separate simulation. 



the spherical surface at the median density, which cor- 
responds to a characteristic drop radius r m . The results 
of numerical calculations in Fig. 2] show both the surface 
tension and surface thickness of the drop tend to con- 
stant values as the drop size increases. For large drops 
the boundary becomes approximately flat so the surface 
thickness becomes independent of drop radius, as does 



the surface tension since its major contribution comes 
from the drop boundary. We now examine the surface 
thickness and tension more carefully in turn. 

To derive an approximate expression for the surface 
thickness we use the analytical approximation Eqn. Q19p 
to the density profile of the inside of the drop, which gives 
the density reduction from the drop center 
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From this, the surface thickness D over which density 
falls from 90% to 10% is given in the large drop limit 
r m > D by 
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(25) 



In the given example in Fig. [3] (12 flavors) this pre- 
dicts that the surface thickness is D m 0.8aQ, which is 
of similar size to the result found by numerical solution 
of Eqn. (fT5)) of ~ Flag, but indicates that the approxi- 
mation for surface thickness in Eqn. (|25|) is not able to 
produce accurate results. The values for surface thick- 
ness of drops found using the many-flavor theory can be 
compared with results from other approximations. For 
the six flavor gas in the large drop limit the many-flavor 
theory approximation Eqn. (|25[) predicts a thickness of 



1.2<2q, and exact numerical integration of the many-flavor 
theory Eqn. (fTS)) predicts thickness 1.6an. The silicon 
six flavor result of Kalia and Vashishta [8[ has a surface 
thickness of 1.6aj5, which is in good agreement with the 
many-flavor result. 

Having used many-flavor theory to predict the density 
profile and surface thickness of an electron-hole drop, 
it is interesting to examine their scaling relationships 
with number of flavors. This is because the scaling re- 
lationships can be experimentally probed [7| by compar- 
ing the surface thickness before and after putting the 
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material under a strain which reduces the valley degen- 
eracy, for example in silicon from six to two flavors. 
These scaling relations will also allow the many-flavor 
results to be further compared with previous theoreti- 
cal work. In Sec. Ill Bl it was shown that the expected 
MFEG uniform density is n ex i/ 8 / 5 , and from Eqn. ([20]) 
Q 2 oc v 2 / b 1 which with Eqn. (f2"5} predicts surface thick- 
ness to scale as D oc v~ x l b . This scaling prediction for 
surface thickness can be compared with numerical results 
for the variation of surface thickness with number of fla- 
vors in Fig. found by solving the differential equation 
Eqn. (|18[) . The coefficient for surface thickness D oc v a is 
a = —0.19995(7) in good agreement with the predicted 
value of —1/5. We can also qualitatively compare our 
scaling result with numerical results [8|, [25J from stud- 
ies of the electron-hole drop in silicon. These studies 
compared results for silicon found at the unstrained six 
flavor with the results at two flavors to attempt to model 
the effect of stress reducing valley degeneracy. Though 
two flavor calculations cannot be accurately given by the 
many-flavor theory, the qualitative variation of surface 
tension and surface thickness should be. The variation 
of surface thickness with number of flavors D ex i/ -1 ' 5 is 
weak, for silicon from six to two flavors the many-flavor 
theory, assuming it is valid, predicts that the thickness 
increases by a factor of 1.2. This compares reasonably 
with the numerical results of Ref. [8j , which predicts an 
increase in surface thickness by a factor of ~ 1.1. 

The dominating contribution to surface energy is at 
the boundary of the drop so the surface tension in large 
drops is approximately the gradient term in the energy 
density, Eqn. (|16p (the main contribution to the surface 
tension) multiplied by the surface thickness D 



(Vn) : 
8n 



D 



it 
8L>' 



(26) 



where we use the additional approximation Vn « n/D. 
Finally, with the relationship found above, D ex ^ _1 / 5 , 
and n oc i/ 8 / 5 found in Sec. Ill Bt this predicts surface 
tension varies with number of flavors as 7 oc v 9 / b . The 
numerical results of Fig. [5) found by solving the differen- 
tial equation Eqn. (fT5|) exactly, predict a coefficient for 
7 oc v a of a — 1.8004(3) in good agreement with the 
analytical result, 9/5. For silicon, reducing the number 
of flavors from six to two, the above result predicts that 



surface tension reduces by a factor of 7. This qualita- 
tively agrees with the variation seen by Refs. [8|, [25| of 
a reduction by a factor of 3, though comparison is dif- 
ficult due to the presence of holes and there being too 
few flavors present for the many-flavor theory to be fully 
applicable. 



V. CONCLUSIONS 

This paper describes a new formalism for calculating 
the behavior of a MFEG. In the many-flavor limit the 
Fermi momentum reduces as pp oc i/ -1 ' 3 so is small com- 
pared with the momenta associated with the strongest 
interactions. Intra-valley interactions are more signifi- 
cant than inter-valley. 

The behavior of a homogeneous MFEG in the 
limit of many-flavors was derived. Specifically the 
exact interacting energy per particle is E[ nt = 
-0.574447(^a* 3/4 m* 1 / 4 ) 71 i/4. ma ki n g it energetically 
favorable for the MFEG to be dense. The formalism was 
found to apply with as few a six flavors over the density 
range 0.03^ < naf < 0.005jA 

The MFEG has short characteristic length-scales 
which motivates a LDA. A gradient expansion of the en- 
ergy density with the lowest order term |Vn| 2 /8n was 
derived, which was applied to electron-hole drops to 
study their density profile and surface properties. Sur- 
face thickness was found to scale as D oc v~ x /^, surface 
tension as 7 oc v 9 / 5 . 

It would be useful to compare our analytical results 
with those from computer simulations to verify our find- 
ings for the uniform MFEG, its polarizability and the gra- 
dient expansion. This would allow the limits over which 
the many-flavor limit applies to be derived more accu- 
rately, and allow the formalism to be applied to more 
physical systems. 
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